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Electroproduction of exclusive 4> vector mesons has been studied with the CLAS detector in the 
kinematical range 1.4 < Q 2 < 3.8 GeV 2 , 0.0 < t' < 3.6 GeV 2 , and 2.0 < W < 3.0 GeV. The scaling 
exponent for the total cross section as 1/(Q 2 +M|) n was determined to be n = 2.49±0.33. The slope 
of the four- momentum transfer t' distribution is = 0.98 ± 0.17 GeV -2 . Under the assumption of 
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sections to be R = 0.86 ± 0.24. A 2-gluon exchange model is able to reproduce the main features of 
the data. 
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I. INTRODUCTION 



Exclusive electroproduction of vector mesons is an es- 
sential tool for exploring the structure of the nucleon and 
the exchange mechanisms governing high-energy scat- 
tering. For low photon virtualities relative to the vec- 
tor meson mass, Q 2 < m v , or in the case of photopro- 
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duction, Q 2 = 0, these processes are well described by 
i-channel exchange of Regge poles (Pomeron, Reggeon) 
— extended objects whose properties can be related to 
the observed hadron spectrum [l|. At high virtualities, 
Q 2 ^> m v , a QCD factorization theorem Q states that 
vector meson production from longitudinally polarized 
photons proceeds by exchange of a small-size system of 
quarks or gluons, whose coupling to the nucleon is de- 
scribed by the generalized parton distributions (GPDs). 
By studying the dependence of exclusive electroproduc- 
tion on Q 2 , one can thus "resolve" the Pomeron and 
Reggeon into their quark and gluon constituents. Ad- 
ditional information comes from the comparison of the 
p°,u> and <fi channels, which couple differently to quarks 
and gluons. The self- analyzing decays of the spin-1 
mesons allow one to study also the helicity structure of 
the 7*iV interaction and, assuming helicity conservation, 
to separate longitudinal and transverse photon polariza- 
tions. 

This article presents data for exclusive (f> vector me- 
son electroproduction off the proton above the resonance 
region, taken with a 5.754 GeV electron beam of the CE- 
BAF accelerator and the CLAS detector at Jefferson Lab 
|3(. The measurement was performed as part of a series 
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of experiments aimed at studying vector meson produc- 
tion in the valence quark region at the highest available 
photon virtualities. The analysis of u> production has 
been completed and the analysis of p production is 
in progress 5]. The analysis of 0- meson production re- 
ported here is based in part on the work of Rcf . Q . 

The 0-meson is unique in that its quark composition is 
mostly ss containing little, if any, u and d flavors which 
populate the valence quarks in the nucleon. Thus, <f> pro- 
duction primarily probes the gluon degrees of freedom in 
the target. High-energy photoproduction of <j> proceeds 
mainly by Pomeron exchange. At large Q 2 , calculations 
based on current GPD models show that the 4> produc- 
tion cross section is dominated by the gluon GPD, with 
only small contributions arising from intrinsic strange 
quarks in the nucleon @,@|. At intermediate Q 2 , a de- 
scription of <f) production based on effective two-gluon 
exchange has been proposed 0], which effectively inter- 
polates between the "soft" and "hard" regimes. Thus, <fr 
production provides us with a clean method of probing 
the gluon field in the nucleon, even at JLab energies. 

A natural framework for discussing exclusive vector 
meson production is the space-time picture in the tar- 
get rest frame (i.e., the laboratory frame) At high 
energies, the interaction of the virtual photon with the 
target proton proceeds by way of fluctuation of the pho- 
ton into virtual hadronic (or quark-antiquark) configu- 
rations that subsequently scatter diffractively off the tar- 
get. This process occurs over a characteristic time given 
by the lifetime of the fluctuation as dictated by the un- 
certainty principle and is given by 
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where v is the photon laboratory energy and Mfl uct is 
the mass of the virtual hadronic state. This interval also 
determines the coherence length in the longitudinal direc- 
tion, Z co h = cAr. In photoproduction or electroproduc- 
tion at Q 2 < m v , this picture is the basis for the success- 
ful vector dominance model (VDM) , where the dominant 
hadronic fluctuations are assumed to be the observed 
ground-state vector mesons (p , </>). Their interaction 
with the target can be described by Pomeron exchange. 
As Q 2 increases, higher-mass states become important. 
Eventually, at Q 2 ^> m v , the fluctuations of the pho- 
ton can appropriately be described as quark-antiquark 
pairs ("dipoles") with transverse momenta k\ ~ Q 2 , or 
transverse size r±_ ~ 1/Q -C I /my. Their interaction 
with the nucleon is described by the gluon GPD, which 
can be interpreted as the "color dipole moment" of the 
target. 

In the context of the space-time picture, measuring 
the Q 2 -dependence of exclusive electroproduction up to 
Q 2 ~ few GeV 2 allows one to vary the transverse size 
of the projectile from "hadronic size" (r± ~ l/my) to 
"small size" (r± ~ 1/Q), thus resolving the structure of 
the target at very different distance scales. At HERA en- 
ergies, where l co h » lfm even for Q 2 ~ few GeV 2 , one 



can neglect the variation of the coherence length with Q 2 
and associate the Q 2 -dependence entirely with a change 
of the transverse size of the projectile. The predictions 
derived in this approximation are nicely confirmed by the 
data, e.g. the decrease of the t-slope with Q 2 , and the 
increase with Q 2 of the exponent governing the energy 
dependence (for a review see Ref. [13). At JLab en- 
ergies, where the coherence length in electroproduction 
is ^coh ^$ 1 fm, we must also take into account its vari- 



ation with Q 



the "shrinkage" of the longitudinal 



size of the virtual photon with increasing Q 2 . Another 
effect modifying the space-time interpretation is the non- 
negligible longitudinal momentum transfer to the target, 
which increases with Q 2 . Nevertheless, the space-time 
picture remains a very useful framework for discussing 
vector meson electroproduction even at JLab energies. 

In the present 0- meson production experiment, the t— 
dependence of the differential cross section was measured 
over a wide range, from the kinematic minimum at t ~ to 
(small CM scattering angle) to t ~ s/2 (large angle). In 
exclusive electroproduction t is related to the transverse 
momentum transfer to the target, Aj_, and thus deter- 
mines the effective impact parameters in the cross sec- 
tion, b± ~ 1/ Aj_ 1 . Exclusive meson production at large 
— (t — to) probes configurations of small transverse size 
in the target. The possibility to vary both Q 2 and t in 
electroproduction allows one to control both the size of 
the projectile and the size of the target configurations 
contributing to the process, and to study their interplay 

Quantitative predictions for the production of vector 
mesons in our kinematic regime have been made by Laget 
and collaborators based on the interactions between con- 
stituent partons (JML model). The high-t behavior of 
the photoproduction cross section of 0-mesons 13| has 
been reproduced using dressed gluon propagators and 
correlated quark wave functions in the proton [14| . Quark 
exchange processes, which contribute also to the pho- 
toproduction of p and lo mesons, have been modeled 
in terms of saturating Regge trajectories. The model 
uses electromagnetic form factors in the Regge amplitude 
[T3,[lH| to describe electroproduction data. However, the 
Q 2 dependence of the 2-gluon amplitude is an intrinsic 
part of its construction, and no additional electromag- 
netic form factors are needed. Therefore, the predicted 
</>-meson electroproduction cross section is parameter free 
and constitutes a strong test of the partonic description 
that underlies the model. The full form for the ampli- 
tudes are given in Refs. [H, H, HH . Thus far, comparisons 
of the JML model for electroproduction have been made 
with lu 0, HH , and p [H, G3] electroproduction data from 



1 The transverse momentum transfer to the target is given by 
= (1 — £ 2 )(i — to), where £ is the fractional longitudinal 
momentum transfer to the target, which in turn is related to the 
Bjorken variable in the kinematics of deep-inelastic scattering, 
i = x B /{2-x B ). 



JLab, and p electroproduction data from HERMES 

One of the leading motivations for the present work is 
the sparse amount of existing <f> electroproduction data. 
The body of (/>-meson electroprodution data at similar 
kinematics consists of early data from Cornell [l9l[20l . [2l| . 
and some data from CLAS at lower energy |22j |. Re- 
cent data on </> electroproduction comes from HERMES 
[H, SI and HERA 0, HE M, H3 at much higher center- 
of-mass energy (W). A summary of the world data in- 
dicating their kinematic range is given in Table [U The 
data from this experiment are complementary to mea- 
surements at collider energies which cover a higher W 
and higher Q 2 range where diffraction mechanisms are 
probed. 

TABLE I: Summary of (f> electroproduction data and kine- 
matic range. 



Experiment 


Q z (GeV 2 ) 


W (GeV) 


Cornell Dixon [19, _20j 
Cornell Cassel [21] 
HERMES [18, 23] 
CLAS [22] 


0.23 - 0.97 
0.80 - 4.00 
0.70 - 5.00 
0.70 - 2.20 


2.9 
2.0 - 3.7 
4.0 - 6.0 
2.0 - 2.6 


HI [24J 
HI [25J 
HI [2JJ 

ZEUS [26] 
ZEUS [26] 


> 7.0 
1.00 - 15.0 
3.00 - 20.0 
7.00 - 25.0 
2.00 - 70.0 


~ 75.0 
40.0 - 130.0 
4.0 - 120.0 
42.0 - 134.0 

35-145 



We have measured 0-meson electroproduction at the 
highest possible Q 2 accessible at CEBAF energies in the 
valence quark regime. The data set covers the kinemat- 
ical regime 1.4 < Q 2 < 3.8 GeV 2 , 0.0 < t' < 3.6 GeV 2 , 
and 2.0 < W < 3.0 GeV. We will present cross sec- 
tions as a function of the momentum transfer —t, the 
azimuthal angle $ between the electron and hadron scat- 
tering planes, as well as the angular decay distributions 
in the rest frame of the 0-meson. Although limitations of 
the statistical sample will preclude determining correla- 
tions between different kinematic variables, the distribu- 
tions will provide insights into the distance scale of the 
interaction and explore kinematics that begin to probe 
partonic degrees of freedom. 



II. KINEMATICS AND NOTATION 



The kinematic variables in exclusive <f> production (see 
Fig.[T]) described by 

e{k)p(P)^e{k')<f>{v)p{P'), (2) 
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Hadron Production Plane (cm.) 



FIG. 1: (Color online) Graphical representation of <^-meson 
electroproduction. Shown from left to right then above are 
the electron scattering plane, the hadron production plane 
and helicity rest frame of the <f) respectively. 

• Q 2 = —q 2 = —(k — k') 2 , the negative four- 
momentum squared of the virtual photon; 

• W 2 = (q + P) 2 , the squared invariant mass of the 
photon-proton system; 

• xb = Q 2 /(2P ■ q), the Bjorken scaling variable; 

• v = P ■ q/M p , energy of the virtual photon; 

• t = (P — P') 2 , the squared four-momentum transfer 
at the proton vertex, is given by 

t = t a -^vX m vtn sin2 (<W2) , where 
in = (E 1 * - E* S 2 - (v 1 * -d* ) 2 

u l u cm ^cmJ \Pcm Pern) 

and the above formulas are calculated using the 
energy and momenta of the virtual photon and (f> 
in the j*p center-of-mass; 

• t' = \t — to\, momentum transfer relative to the 
kinematic limit —to, which increases with Q 2 and 
decreases with increasing W; 

• The coordinate system in the j*p center-of-mass 
is defined with the z-axis along the direction of 
the virtual photon, and the y-axis normal to the 
hadronic production plane along pj m x p f m ; 

• $, the angle between the hadron production ("i*4>p) 
plane and the electron scattering (ee'7*) plane fol- 
lowing the convention in Ref. [28j| 2 : 



are k, k' , P, P' and v which are, respectively, the four- 
momenta of the incident electron, scattered electron, tar- 
get proton, scattered proton and the (/>-meson: 



2 The azimuthal angle <E> used here is — <j> from the "Trento con- 
vention" |29( . 
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• cos Oh and Sh, decay angles of the K + in the he- 
licity frame [28|, which is defined in the rest frame 
of the </>-meson with the z-axis along the direction 
of the 0-meson in the 7*p center-of-mass system; 

m tp = (f>jj — <f>, azimuthal angle that simplifies the 
angular decay distributions when s-channel helicity 
is conserved (SCHC). 

The electroproduction reaction integrated over the de- 
cay angles of the 0-meson can be described by the follow- 
ing set of four independent variables: Q 2 , —t, $ and W. 
For the analysis of the decay distribution, the additional 
variables cos Oh and i\) are required. In total there are six 
independent variables in the approximation of negligible 
<j> width. 

III. EXPERIMENT 

The experiment was conducted with the CEBAF Large 
Acceptance Spectrometer (CLAS) Q located in Hall B 
of the Thomas Jefferson National Accelerator Facility. 
The CLAS spectrometer is built around six independent 
superconducting coils that generate a toroidal magnetic 
field azimuthally around the beam direction. The az- 
imuthal coverage is limited by the magnetic coils and 
is approximately 90% at large angles and narrows to 
50% at forward angles. Each sector is equipped with 
three regions of multi-wire drift chambers and time-of- 
flight counters that cover the angular range from 8° to 
143°. Charged-particle trajectories are tracked through 
the field with the drift chambers, and the scintillators 
provide a precise determination of the particle flight time. 
In the forward region (8° to 45°), each sector is fur- 
thermore equipped with gas-filled threshold Cerenkov 
counters (CC) and electromagnetic calorimeters (EC). 
The Cerenkov counters are used to discriminate electrons 
from pions, and the calorimeters are used to measure the 
energy of electrons and photons. 

The data were collected between October 2001 and 
January 2002 with a 5.754 GeV electron beam incident 
on a 5 cm- long liquid hydrogen target. The typical beam 
current was 7 nA. The CLAS torus magnet was set to 
3375 A with a polarity that caused negatively charged 
particles to bend in towards the beamline. The inclusive 
electron trigger fired when signals in the forward elec- 
tromagnetic calorimeter exceeded a predefined threshold 
in coincidence with a hit in the Cerenkov counters. The 
kinematical domain of the selected sample corresponds 
approximately to Q 2 from 1.5 to 5.5 GeV 2 and W be- 
tween 2 and 3 GeV. The typical experimental dead time 
was about 8% with a trigger rate of about 1.5 kHz. 

IV. EVENT RECONSTRUCTION 

The (/)-mesons were detected using the charged-particle 
decay mode into K + and a K~ . Events corresponding 




FIG. 2: Energy deposited by the electron candidates in the 
electromagnetic calorimeter versus momentum. The lines 
show the selection cuts for good electrons as described in the 
text. 

to ep — > epK + (K~) were classified initially by requiring 
at least one negative track and two positive tracks. Nor- 
mally the K~ remained undetected due to the limited 
acceptance for negative particles at this high magnetic 
field setting. After calibration of the spectrometer, the 
momentum of each particle was determined with a frac- 
tional resolution of about a percent using the track seg- 
ments in the drift chambers. The momentum resolution 
is sufficient to identify the missing particle as a K~ . 

The identification of good electrons is the crucial first 
step and is accomplished through energy and momentum 
cuts |6j . After selection of tracks within the fiducial vol- 
ume of the detector, the momentum of the electron can- 
didate track in each event was required to correspond to 
the energy deposition in the electromagnetic calorimeter 
and the visible energy be greater than 0.2 GeV (Fig. [2]). 
Pions were rejected by requiring a minimum energy of 
0.06 GeV in the inner layer of the calorimeter and a pulse 
height in the Cerenkov counter corresponding to at least 
2.5 photoelectrons [30l l3ll| . 

The two positive tracks in the fiducial volume were 
identified as a proton and K + using the measured flight 
time {5T ~ 160 ps) from the target to the time-of-flight 
counters [13], a typical distance of about 5 m. Fiducial 
volume cuts were made to cut out tracks in inefficient 
parts of the detector and small momentum corrections 
were applied to compensate for uncertainties in the mag- 
netic field and detector positioning. The time of the in- 
teraction was determined using the vertex time of the 
electron corrected to the time of the bunch crossing of 
the machine. Using the known momenta of each of the 
tracks, the vertex time was computed making assump- 
tions for the mass of the particle and comparing to the 
time of the bunch crossing. Events were kept where the 
two positive tracks were consistent with the assignment 
of one proton and one K + . Tracks were identified as pro- 
tons when the projected vertex time assuming a proton 
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FIG. 3: Mass computed from the flight time versus momen- 
tum for positive particles. The top band corresponds to pro- 
tons, the middle band corresponds to K +, s and the lower 
enhancement at 1.5 GeV/c momentum is due to pion con- 
tamination. 



mass differed from the interaction time by less than 0.75 
ns, and as a K + when the projected time assuming a 
kaon mass differed from the interaction time by less than 
0.6 ns. In cases where one track satisfied both criteria, 
the ambiguity was resolved using the second track. The 
number of events where both tracks satisfied both crite- 
ria was less than 1% and were eliminated. The calculated 
mass versus momentum, shown in Fig. 02 indicates that 
at high momenta, there remain a number of pions that 
are identified as kaons in the sample. 

Once the electron, proton and K + tracks were iden- 
tified, the missing mass was computed and is plotted in 
Fig. 2] A clear peak is found at the mass of the K~ , which 
corresponds to the exclusive reaction ep — > epK + K~ . A 
2cr cut was applied to the epK + X events to select the 
sample of interest. For those events, the four-vector for 
the K~ was constructed by setting the three-momentum 
equal to the missing momentum of the ep — > epK + X re- 
action, and the energy was then calculated usin g th e K~ 
mass recommended by the Particle Data Group [33j . The 
fraction of events where the K~ was detected in the de- 
tector was so small that they were not treated differently 
than the rest of the sample. 



A. (f> Event Identification 

The sample satisfying the epK + (K~) criteria contains 
27,950 events out of 947,300 epK+X candidates. The 
sample includes all physical processes that contribute to 
this final state, as well as real fts and background from 
misidentified pions. Fig. [5] shows the K + K~ invariant 
mass (Mkk) for the entire data set with a clear 0- meson 
peak. This distribution is simultaneously fit to a Gaus- 
sian plus an empirical phase space function for the back- 



FIG. 4: (Color online) Distribution of epK + X missing mass. 
The vertical lines indicate the cuts placed to select events with 
a missing K~ . 
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FIG. 5: (Color online) K + K invariant mass including all 
data cuts and a fit to <f> peak with Eq.[3] 



ground, 

FIT = AG(a, ft) + B 1 yjMu?-M& + B 2 (lV1 kk 2 - M 2 h ) , 

(3) 

where G(cr, //) is a Gaussian distribution, M t ^=0.986 GeV 
is the threshold for two kaon production and A, B±, and 
Bi are parameters of the fit. This fit yields = 792±52, 
a mean n = 1.0194 ± 0.0005 GeV, and a width of 
a — 6.5 ± 0.6 MeV. The signal-to-background ratio for 
this fit was 0.56. The mean and width are fixed to these 
values for all subsequent fits to the invariant mass dis- 
tribution to constrain the fits with limited statistics in 
specific kinematic bins. A total of 37 distributions were 
fitted to extract the <fi signal in various kinematic bins 
(see subsequent sections for details) . The average x 2 per 
degree of freedom for all the fits was 1.07, indicating that 
deviations from the fit function are statistical in nature. 
There are competing physics channels that also lead to 
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A(1520) 





y 



K + 




A*(L*) 
P P 



FIG. 6: (Color online) Sketch of exclusive 0(1020) production 
via Pomeron exchange and of excited hyperon production, of 
which A(1520) is an example. This is the primary physics 
background for 0(1020) production. 



the same final state. The majority of these backgrounds 
come from the production and subsequent decay of high- 
mass hyperons produced via ep — ► e A" + A*(£*) as illus- 
trated in Fig. [6] The Dalitz plot in Fig. [7] clearly shows 
the dominant A(1520) background contribution (horizon- 
tal strip), as well as the 0(1020) (vertical strip). There 
are additional contributions from the higher-mass states 
such as A*(1600), A* (1800), A* (1820), E*(1660), and 
£*(1750) but they cannot be separately identified. In or- 
der to avoid the introduction of holes in the acceptance, 
no cuts are made to remove these hyperon backgrounds. 
Instead they are taken into account during the fitting 
procedure by assuming they contribute to the smooth 
background under the 0-meson peak. Nevertheless, many 
different fits were performed removing events in the peak 
of the A* (1520) to study this systematic with no indica- 
tion that they changed the results significantly. These 
studies focused on the t-distributions, since the effective 
momentum transfer in A* reactions is very flat compared 
to that expected from 0-meson production. 



ACCEPTANCE CORRECTIONS 



Particle interactions and event reconstruction in the 
detector were simulated using a GEANT-based Monte 
Carlo called GSIM [34]. The events were gen erated ac- 
cording to a VDM-inspired cross section [l0| with the 




1.8 2 

M 2 K+K -(GeV 2 ) 



FIG. 7: Dalitz plot of M% K versus 

Mkk- The well-defined 
horizontal strip is the A(1520) band. The vertical strip is the 
0(1020) band. 



following form: 



T {Q\W) 



a^W) [1+Re] 
(l + Q 2 /M 2 f 

(W 2 



M 2 ) exp( 



-bf) 



^{W 2 ~M P 



=(4) 



AE e {E e -v)-Q 



2)2 +4M /2Q2 
2 



2v 2 + Q 2 



(5) 



where <r$(0, W) is the (transverse) photoproduction cross 
section, E e is the incident electron beam energy, e is the 
virtual photon polarization parameter and R is the ra- 
tio of the longitudinal to transverse cross section. The 
parameters of the model were tuned during preliminary 
analysis and found to reproduce the general features of 
the data. The main variation from the conventional 
VDM model was in the propagator, where preliminary 
data seemed to indicate a stronger dependence on Q 2 
and an exponent of 3 was used instead of 2. 

The acceptance function is a combination of the geo- 
metrical acceptance of CLAS, the detector efficiencies of 
the scintillators and drift chambers, the track reconstruc- 
tion efficiency, and the event selection efficiency. The 
Cerenkov detector [30( is not well modeled in GSIM, and 
its efficiency was determined separately using the data. 

The acceptance was defined in each bin of a 6- 
dimensional table as the ratio of reconstructed to gener- 
ated Monte Carlo events. In order to account for corre- 
lations between all kinematic variables, a total of 33,600 
acceptance bins are defined in the kinematic variables 
Q 2 , —t, W, $, cos6h and 0. The binning selection is 
given in TablclTTl for the first three variables and uniform 
binning was used for $ (6 bins), cos 6h (5 bins) and ip (8 
bins). The projected 2-D acceptance surface in Q 2 and 
—t and the 1-D projections in Q 2 , t, and W are shown in 
Fig. [8] . The projected 2-D acceptance surface in cos 6h 
and i/j is shown in Fig.[9j as well as the 1-D projections 
in cos9h, 0, and The variation of the acceptance is 
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FIG. 8: (Color online) 2-D Acceptance in Q 2 and t, as well 
as the 1-D acceptance in Q 2 , W, and —t. Error bars are not 
shown; the lines are present to guide the eye. The axes in the 
2-D plot in a) have the same range as that of the axis of the 
two 1-D plots in b) and d). 




FIG. 9: (Color online) 2-D Acceptance in cos#h and ip, as 
well as the 1-D acceptance in cos Oh, and ip. Error bars 
are not shown; the lines are present to guide the eye. The 
axes in the 2-D plot in a) have the same range as that of the 
axis of the two 1-D plots in b) and c). 



relatively smooth as a function of these variables (except 
for $, which is a reflection of the CLAS torus coils) and 
is of the order of 1-3%. 

Events that fell into bins with extremely small accep- 
tances (< 0.2%) were eliminated to avoid biases due to 
statistical fluctuations in those bins. The losses were es- 
timated and corrected by using the ratio of Monte Carlo 
acceptance-weighted events to generated events. 





No. 


Bin Definition 


Q 2 


5 


1.4-1.8 


1.8-2.2 


2.2-2.6 


2.6-3.0 


3.0-3.8 






w 


4 


1.9-2.1 


2.1-2.5 


2.5-2.7 


2.7-2.9 








-t 


6 


0.4-0.8 


0.8-1.2 


1.2-1.6 


1.6-2.0 


2.0-2.4 


2.4-3.6 




t> 


7 


0.0-0.2 


0.2-0.4 


0.4-0.6 


0.6-0.8 


0.8-1.0 


1.0-2.0 


2.0-3.6 



TABLE II: Binning for the acceptance calculation in Q 2 
(GeV 2 ), -t (GeV 2 ), and VF(GeV). An additional acceptance 
table was also generated for t' in the place of t, but it is not 
an independent variable. 



VI. 



RADIATIVE CORRECTIONS 



The radiative effects were calculated in two distinct 
steps. The external radiative process, which is the finite 
probability that the incoming or scattered electron will 
radiate a hard photon in the presence of a nucleon in the 
target other than the one associated with the event, is 
taken into account during the Monte Carlo acceptance 
calculation. The internal radiative corrections include 
the Bremsstrahlung process for the incoming or scat- 
tered electron in the presence of the nucleon associated 
with the event, as well as diagrams such as vacuum po- 
larization, which are not accounted for during the ac- 



ceptance calculation. These are included in the correc- 
tion factor F rac i using the radiative correction code EX- 
CLURAD setting the controlling parameter v cu t = 0.047 
GeV 2 [35J . F ra d is calculated in each W and <I> bin as the 
ratio (Jrad/ v-norad (variable 6 in Eq. 75 from Ref. [35|). 
where a nora d is the cross section calculated without any 
radiative effects (i.e. the Born cross section) and a ra d 
is the cross section calculated with radiative effects in- 
cluded. The correction factor for various W bins is shown 
as a function of <& in Fig. [10] The correction was com- 
puted in bins of W and «3> for average values of Q 2 and 
cos Oqm (directly related to —t) because the correction 
was found to change less than 2% over the range of Q 2 
and cos 9cm [III . 



VII. CROSS SECTIONS 



The reduced 7*p — 
is given by 

<y(Q\w) 

T{Q 2 ,W,E e ) 



electroproduction cross section 



da 



T{Q 2 ,W 1 E e )dQ 2 dW 



W(W 2 



M 2 ) 



4tt M 2 E 2 Q 2 1-e' 



(6) 



where T(Q 2 , W, E e ) is the virtual photon flux factor. We 
can extract the <h cross section from the data via 



da 



1 



n w 



dQ 2 dW Br £ mt AQ 2 AW' 



(7) 



where Br is the branching fraction ratio of <p — > 
K+K- = 0.491 ± 0.009 [H, £ int =2.98 x 10 40 cm" 2 is 
the live-time-corrected integrated luminosity, AQ 2 and 
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100 200 300 

4> (degrees) 

FIG. 10: (Color online) Plot of radiative correction F ra d as a 
function of $ for assorted values of W from 2.0 to 3.0 GeV. 
The correction for each W value was computed for < Q 2 >= 
2.47 GeV 2 and < cos# c . m . >= 0.345. 



Q 2 range 


W range 


1.4 < Q 2 < 1.8 


2.10 < W < 2.90 


1.8 < Q 2 < 2.2 


2.10 < W < 2.90 


2.2 < Q 2 < 2.6 


2.10 < W < 2.90 


2.6 < Q 2 < 3.0 


2.10 < W < 2.70 


3.0 < Q 2 < 3.8 


2.10 < W < 2.70 



TABLE III: W range for each Q 2 bin. 



cross section was fit to the function 



to determine the scaling behavior. For this data we de- 
termined the parameter n — 1.97 ± 0.84. The measured 
exponent spans the range expected for the dependence 
on Q 2 due to VDM (n — 2) to hard scattering (n = 3 for 
fixed momentum transfer t). 



AW are the corresponding bin widths modified appro- 
priately when not completely filled due to kinematics, 
and n\Y is the result of a fit to the Mkk distribution 
weighted by acceptance, CC efficiency correction and ra- 
diative corrections. The binning in Q 2 , —t, t' , and W 
for the extraction of the cross section in those variables 
is shown in TableHTl We emphasize here that we have 
performed a fit to Eq.[3] to determine the signal nw and 
the estimated background under the peak for each entry 
in the table. The differential cross section in a variable 
X is given as 



da 
dX 



AX 



(8) 



where a(Q 2 ) is the total cross section in a bin in X. The 
cross sections presented in this paper have been corrected 
for the bin size and are quoted at the center of each bin. 



A. Total Cross Section a(Q 2 ) 

The cross section a(Q 2 ) as a function of Q 2 is ob- 
tained by integrating over W due to the limited statistics. 
Each event was weighted for acceptance, radiative effects, 
the CC efficiency, as well as the virtual photon flux fac- 
tor. The invariant mass distribution (Mkk) of weighted 
events in each Q 2 bin was then fit to Eq.[3] The bins 
used in the analysis are given in Tablc lTTTl The range in 
W was restricted at the low end where acceptance correc- 
tions change rapidly and are large, and at the high end to 
match the high end of the kinematically accessible range. 
The cross section for each of the bins was calculated ac- 
cording to Eqs.[6]and[7] A small correction (~l-2%) was 
then applied to adjust the bin-averaged cross section to 
the center of the bin (37[. The values for the cross sec- 
tion in each Q 2 bin are shown in Table ftVl The total 



Q 2 (GeV 2 ) 


< e > 


a (nb) 


1.6 


0.488 


9.9 ± 3.2 


2.0 


0.479 


10.4 ± 2.5 


2.4 


0.471 


6.7 ± 1.9 


2.8 


0.464 


5.9 ± 2.4 


3.4 


0.452 


3.6 ± 1.8 



TABLE IV: Total cross section a(Q 2 ) and kinematics of each 
data point, along with the center of each Q 2 bin. < e > is the 
average virtual photon polarization in each bin. 



B. Differential Cross Section in t', da/dt' 

The differential cross section in t' was extracted in 
seven bins in £' by fitting Eq.[3]to the K + K~ mass distri- 
bution to determine the 4> signal and background in that 
particular bin. The average x 2 per degree of freedom for 
these fits was 1.2. The signal-to-background ratio varied 
from bin to bin, ranging from 0.33 to 0.86. The lowest 
signal-to-background ratio occurred in the mid range of 
t'. The resulting values for the cross section in each t' 
bin are shown in TablcfVl 

In cases of limited statistics, da/dt' is often used in- 
stead of da/dt in order to eliminate kinematic corrections 
due to — to, which varies with Q 2 and W. This procedure 
is most useful when the cross section factorizes into terms 
that depend only on t and terms that depend on Q 2 and 
W, aside from the threshold dependence, as in the VDM 
model. Indeed, our measured differential cross section in 
t' show very similar trends as previous data, namely they 
are consistent with diffractive production (e -6 *'* ') [23 |. 
Fig. HU shows an exponential fit to the measured differen- 
tial cross section, which yields a = 0.98 ±0.17 GeV -2 . 
At high energies, the slope can be directly interpreted in 
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0.5 1 1.5 2 2.5 

t' (GeV 2 ) 

FIG. 11: (Color online) Plot of da/dt' along with an expo- 
nential fit. 



reliably, so that bin was dropped. The first bin reported 
contains a significant correction, but was included with 
an increased systematic error. Subsequent bins had small 
or no corrections. The values for the cross section in each 
—t bin are given in Table fVTl 



-t (GeV") 


da/dt (nb/GeV*) 


0.6 


10.7 ± 3.1 


1.0 


0.8 ± 1.0 


1.4 


3.4 ± 1.0 


1.8 


1.0 ± 0.5 


2.2 


1.4 ± 0.5 


3.0 


0.5 ± 0.2 



TABLE VI: Differential cross section da/dt and kinematics of 
each data point, —t is the center of each bin at an average 
value of e = 0.47. 



terms of the transverse size of the interacting configura- 
tion, as described later when presenting results. In that 
limit, the small value of the exponential slope implies the 
interaction takes place at very short distances inside the 
nucleon. 



H (GeV 2 ) 


da/dt' (nb/GeV") 


0.1 


9.4 ± 2.9 


0.3 


4.4 ± 1.9 


0.5 


3.1 ± 1.2 


0.7 


2.7 ± 1.2 


0.9 


4.0 ± 1.3 


1.5 


1.6 ± 0.4 


2.8 


0.5 ± 0.2 



TABLE V: Differential cross section da/dt! and kinematics of 
each data point, t! is the center of the bin, and corresponds 
to an average value of e = 0.47. 



C. Differential Cross Section in t, da/dt 

The differential cross section is easiest to compare with 
theory if it is computed in terms of the Mandelstam vari- 
able t. The cross section is given as 



da 
~dt 



<J{Q 2 



At ■ Corr(^o)' 



(10) 



where At is the bin size and Corr(to) is a correction factor 
to account for the fact that the kinematic limit to(Q 2 , W) 
varies across the bin. The yield was extracted over the 
ranges of Q 2 and W given in Table HIIl in six bins in —t. 
The kinematic threshold to varies between -0.09 and - 
1.14 GeV 2 for extreme values of Q 2 and W. For the 
bin corresponding to < — t < 0.4 GeV 2 , the threshold 
varies so much that corrections could not be modeled 



D. Differential Cross Section da/d& and test of 
SCHC 

The cross section dependence on the angle $ between 
the electron and hadron scattering planes takes the fol- 
lowing form: 



da 



a + ecTTT cos 23> + \/2e(l + e)cLT cos $ 



(11) 



where <jlt and <jtt are the interference terms between 
the longitudinal and transverse contributions to the cross 
section. If helicity is conserved in the s-channel (SCHC), 
then both of these terms will vanish. The magnitude of 
these interference terms can therefore be used as a test 
for the validity of SCHC. 

The differential cross sections in $ were extracted in 
the same manner as the other differential cross sections 
(Eq.[5]) after integrating over Q 2 , —t and W. The cross 
section da/dQ was extracted in six bins in <!>. The cross 
sections, along with a fit to Eq.QTJ are shown in Fig.[L2"l 
The fit yields a value of <ttt = —1.1 ± 3.1 nb and &lt = 
2.2 ± 1.1 nb with a chi-squared per degree-of-freedom of 
X 2 /D.F. = 1.3. A fit of the d<j/d<$> distribution to a 
constant, constraining the interference terms to be zero, 
yields a x 2 /D.F. = 1.6. The small change in the goodness 
of fit between the two cases leads us to conclude that the 
precision of this experiment is insensitive to violations of 
SCHC for </>-meson production in our kinematic domain. 



VIII. ANGULAR DECAY DISTRIBUTIONS 

The angular decay distribution of the K + in the (j> 
rest frame describes the polarization properties of the 
</>-meson. The scattering amplitude for vector meson 
electroproduction 7* + N — > P + V can be expressed in 
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FIG. 12: (Color online) der/d<£> vs The curve shows a fit to 
Eq.[TTJ which is used to determine gtt and olt- The dotted 
line is a fit to a constant function which is expected from 
SCHC. 



terms of the helicity amplitudes T\ v \ p \ ,\ N , where is 
the helicity of each particle (i=V, P, 7, N). The vector 
meson spin density matrix is derived from these helicity 
amplitudes by exploiting the von Neumann formula 



1 



p{V) = 7r7>( 7 *)Tt, 



(12) 



where p("f*) is the spin-density matrix of the virtual 
photon. The details of this derivation can be found in 
Ref. [HI]. The density matrix element is denoted p"-, 
where the index a can be related to the virtual photon 
polarization, a = — 2 for purely transverse photons, 
a = 4 for purely longitudinal photons, while other values 
correspond to longitudinal-transverse interference terms. 
The indices ij correspond to the helicity state of the vec- 
tor meson [38] . In cases where the data do not allow for 
a ul /&t separation, the unseparated matrix elements rfj 
can be parameterized as: 



„04 



Pi 



eRpfj 



eR 



rg = VRj 



1 + eR 



a 



0-3 



eR 



a = 5 



(13) 
(14) 
(15) 



Recall that R is the ratio of longitudinal to transverse 
cross section. The angular distribution of the K + is usu- 
ally described in the helicity frame, defined in the rest 
frame of the 0-meson with the z-axis oriented along the 
<^-meson in the 7*p center-of-mass. The full decay dis- 
tribution, which we denote by Wf(cos0jj, $, 4>h ), can be 
found in the literature (39j , but will only be given here in 
simplified forms. In particular, further analysis of angu- 
lar distributions is done under the assumption of SCHC, 
which leads to considerable simplifications with the in- 



troduction of ip = 4>h — <& and the following constraints 

a/Rcos<5 



-Imr 
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' 1-1 



= Rer°i = 
-lmr 2 i^i = 

'00 — 



V8(l + ei?V 
1 

2{l + eR) ] 

eR 
1 + eR' 



(16) 
(17) 
(18) 



All other rg-'s are 0, and v Re is the ratio of the longitu- 
dinal to transverse amplitudes. The angular distribution 
becomes a function of two variables only and is given by: 



W (cob 9 H ,ip) 



1 



8?r 1 + eR 
2(1 + eR)e(r{_!) sin 2 9 H cos 2^ 



sin 2 H + 2ei?cos 2 H 



+ A(l + eR) y/e(l + e)(Rer 5 10 )sm28 H cos^ .(19) 

In order to extract the rfj parameters from the mea- 
sured angular distribution, we use two 1-dimensional pro- 
jections of the full angular distribution. 



A. Polar Angular distribution projection 

To obtain the polar angular distribution, an integra- 
tion of the full angular distribution Wp over <j)H yields 



W (cos H ) 



(l-r ^) + (3r 4 -l)c<^^ ; 



(20) 



which is independent of SCHC. In order to obtain this 
projection from the data, the K + K ~ invariant mass dis- 
tribution is plotted in five bins in cos Oh (0.40 units of 
cos Oh each) . The same fit to a Gaussian plus a polyno- 
mial background was made to extract the weighted yields 
in each of these bins. The fit to da/dcosOn in Fig. [T3l 
yields a value rj$ = 0.33 ± 0.12 with a x 2 /D.F. = 1.7. 
With the additional assumption of SCHC, this parame- 
ter can be used to determine the ratio of the longitudinal 
to transverse cross sections as 



R 



,,04 
'00 



1.05 ±0.38, 



(21) 



e(l - rgg) 

where we have used the average value of < e >= 0.47. 

B. Angular distribution projection in ip 

After an integration of Wp in cos 6^, a substitution 
of <j>H = t/j + and an integration in <!>, the projected 
angular distribution in tp is given as 



W{il>) 



1 

2^ 



l + 2e{r\_ 1 ) cos 2ip 



(22) 
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FIG. 13: (Color online) Unnormalized polar angular decay 
distribution of the K + integrated over all Q 2 values plus a fit 
to Eg, 1201 Also shown is the extracted rgo parameter. 
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FIG. 14: (Color online) Unnormalized azimuthal angular dis- 
tribution extracted for all Q 2 values plus a fit to Eg. 1221 The 
value of r\_ 1 can also be used to determine R. 



of backgrounds and their uncertainties were also limited 
by the finite sample size. The total systematic error of 
18.6% was added in quadrature with the statistical errors 
in all quoted cross sections. 

The acceptance correction contributes to the system- 
atic error in two distinct ways. The uncertainty of 6% 
introduced by eliminating events with very large weights 
(i.e. very low acceptance) was estimated by changing the 
maximum weight allowed and recomputing the extracted 
cross section. The uncertainties introduced by the use of 
our acceptance table (12%) were estimated by combining 
bins and comparing the extracted result to the average 
of the constituent bins. 

To estimate the systematic uncertainty due to the un- 
known distribution of backgrounds, the functional form 
of the background (see Eq.[3]) was modified by adding a 
term proportional to {M^ K — M^ h ) and refitting the —t 
and Q 2 distributions. The new fits were less constrained, 
but the average change in cross section was 9%. We found 
that the extraction of the slope parameter was fairly 
robust to these changes. In addition, the fitted invari- 
ant mass distributions included some background due to 
misidentified pions. The estimated uncertainty due to 
this contamination under the peak was estimated to be 
7%. 

The systematic uncertainty in the placement of the cut 
to select the K~ from the epK + missing mass (5%) was 
investigated by varying the cut and observing the effect 
on the cross sections. The systematic error associated 
with the bin centering correction is almost negligible (~ 
1%). The contribution to the systematic error from the 
radiative correction was estimated to be ~ 3 % and is 
described in more detail in Ref. [35|. The fluctuation 
in the number of photoelectrons in the Cerenkov counter 
over the course of the run can cause a systematic error in 
the Cerenkov counter efficiency correction. This leads to 
a systematic of < 1%. Finally, the procedure to estimate 
the correction due to the to kinematic cutoff in the first 
i-bin introduces a 25% systematic error in that bin. 



which assumes SCHC. The factor of 1/2tt is a normal- 
ization factor. A fit of da /dip to Eq.[22] is shown in 
Fig.[H The fit yields a value = 0.38 ± 0.23 with a 
X 2 /D.F. = 1.3. The ratio of longitudinal to transverse 
cross sections can also be computed from r\_ 1 (Eq.[T7|) 
and gives R = 0.72 ± 0.3, in agreement with the value 
obtained previously. 



IX. SYSTEMATIC UNCERTAINTIES 



Source 


Aa % 


Acceptance correction 


13.4 


Background functional form 


9.1 


Misidentified pion background 


7.0 


epK+(X) cut 


5.0 


Bin-centering correction 


1.0 


Radiative correction 


3.0 


Cerenkov efficiency correction 


1.0 


Total 


18.6 



TABLE VII: Table of systematic errors. 



The relatively low number of measured eft events causes 
statistical errors to dominate. The sources of systematic 
errors in this experiment are summarized in Table lViTl 
The major sources of systematic errors are due to accep- 
tance corrections and estimation of backgrounds. Studies 



X. DISCUSSION 

The measurements of <r(Q 2 ) from the present analysis 
are shown along with other data on (ft electroproduction 
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FIG. 15: (Color online) Total cross sections as a function of 
Q 2 for our data (red full circles), previous JLab data (open 
circle s) 1221 ] . Cornell data (stars) for W between 2 and 3.7 
GeV &~HERMES data (triangles) for W between 4 and 
6 GeVH, and HERA data (squares) at high W 0. The 
curves show the predictions of the JML model at W=2.9, 2.45 
and 2.1 GeV (top to bottom). 



[H [U, HI [H HI in Fig.HH The one overlap point at 
Q 2 = 1.5 GeV 2 is in good agreement with the previous 
CLAS measurement (22J. The data sets span the range 
from threshold at W=2 GeV up to HERA energies. 

The data sets have a similar trend as a function of 
Q 2 and increase monotonically as a function of W. The 
three curves using the JML model at W = 2.1, 2.45 and 
2.9 GeV are also plotted for Q 2 greater than 1.5 GeV 2 . 
The calculation for W—2A5 GeV, which is close to the 
average of our data, seems to overestimate our data by 
about a factor of two, although it does reproduce the 
existing Cornell data from Ref. [2lJ. The Cornell data 
has a much wider acceptance range in W between 2.0 
and 3.7 GeV, so in fact it could be representative of the 
cross section at higher W. The new data from CLAS, 
together with the existing world data, in particular the 
data from HERA, indicate that the qualitative behavior 
as a function of Q 2 does not change between threshold 
and a W of about 100 GeV. 

Of interest is the applicability of factorization and the 
formalism of GPDs to meson production in general, and 
</> production in particular. QCD factorization makes 
certain asymptotic predictions about the cross section, 
namely that the longitudinal part of the cross section, 
ctl, becomes dominant as Q 2 increases, and that the dif- 
ferential cross section will scale as 1/(Q 2 ) 3 at fixed t and 
xb- For a slow variation of the cross section over the 
range of xb of the data (0.2-0.5), this prediction can be 
compared to the Q 2 dependence integrated over W and 
t, although quantitative estimates are modified by power 
corrections as well as kinematics near threshold. On the 
other hand, the VDM model predicts the cross section to 
scale as 1/(Q 2 + M 2 ) n with n = 2. The Q 2 range of our 




Q 2 (GeV 2 ) 

FIG. 16: (Color online) Fit to the cross section as a function 
of Q 2 distribution to determine scaling using data from the 
present experiment and CLAS data from Ref. 



data is limited, but in combination with previous CLAS 
data at lower Q 2 [HJ (see Fig.QU) we can determine the 
scaling exponent of 1/(Q 2 + M 2 ) n to be n = 2.49 ±0.33. 

Present theoretical calculations of the cj> production 
cross section based on GPD models suffer from consid- 
erable quantitative uncertainties when applied to fixed- 
target energies. At HERA energies the approach taken 
in Ref. [iof . which relies on the equivalence of leading- 
order QCD factorization with the dipole picture of high- 
energy scattering, gives a good description of the abso- 
lute cross section, as well as of subtle features such as the 
change of the W— and i-dependence with Q 2 . Essential 
for the success of this approach is the fact that the effec- 
tive scale of the gluon GPD, Q 2 $, is considerably smaller 
than the external photon virtuality, Q 2 , as has been con- 
firmed by detailed quantitative studies [4l[. The same 
is expected in vector meson production at fixed-target 
energies; however, implementing it in a consistent man- 
ner in these kinematics has so far proven to be difficult. 
Leading- twist, leading-order QCD calculations of the (f> 
production cross section at JLab and HERMES energies 
done with the assumption that Q 2 S = Q 2 [7| overestimate 
the measured cross section by a factor 5-10 and predict 
too steep an energy dependence. A satisfactory solution 
to this problem likely requires a comprehensive approach 
that combines contributions from small-size (~ 1/Q) 
and hadronic-size configurations in the virtual photon 
at moderate coherence lengths (ct < lfm), and possi- 
bly higher-order (NLO) QCD corrections. We note that 
a modified perturbative approach Q which includes the 
intrinsic transverse momentum in the meson wave func- 
tion has been fairly successful in reproducing the mea- 
sured cross sections down to relatively low Q 2 and W . 

The four-momentum transfer distribution probes the 
size of the interaction volume. At high energies, the 
exponential slope (see Fig fTTj) is directly related to the 
transverse size ~ \R 2 nt in analogy to the classical 
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FIG. 17: (Color online) Exponential slope 6^ plotted as a 
function of the fluctuation parameters cAr for the world data. 
The data at high W measure an asymptotic slope correspond- 
ing to long fluctuation times. At low W and relatively large 
Q 2 , the fluctuation times becomes small and constrain the 
size of the interaction volume. 



scattering of light through an aperture of radius Ri n t ~ 
0.38 fm. At energies close to threshold, as in the present 
experiment, this interpretation requires some modifica- 
tion. When the coherence length becomes comparable to 
the size of the target, longitudinal shrinkage occurs and 
this also causes a decrease of the exponential slope. The 
longitudinal size is related to the fluctuation time At of 
the virtual meson, which can be estimated through un- 
certainty principle arguments, and is given by Eq.[T] The 
nature of the interaction becomes more point-like as Q 2 
increases and the fluctuation time decreases. This tran- 
sition should be observed as a decrease in the measured 
slope parameter. Since the differential cross section in t' 
was extracted for all Q 2 , the value for corresponds to 
the average value of cAr=0.46 fm. The slope parame- 
ters for various experiments are shown in Fig.[T7| for the 
world data on <f> electroproduction. The measured slower 
fall-off of the ^distribution, corresponding to the small 
slope parameter, is consistent with the expectation that 
short interaction time probe small ss dipoles. 

The differential cross section in —t is compared to the 
JML model in Fig.QU The data covers 1.4 < Q 2 < 3.8 
GeV 2 and the JML model predictions [12] are plotted for 
fixed values of Q 2 from 1.6 to 5 GeV 2 . The data tend to 
have a shallower slope than the calculation, but there is 
general agreement. This agreement is highly non-trivial 
since the few parameters of the model have been fixed 
at the real photon point and kept frozen in the virtual 
photon sector. Our data confirm both the Q 2 and —t 
dependence of the cross section that are built into the 
dynamics of the ss loop and the 2-gluon loop. 

The angular decay distributions provide information 
on the longitudinal part of the production cross section. 
We have extracted values of cttt and &lt from the cross 



FIG. 18: (Color online) da/dt vs —t for the entire Q 2 range 
and the JML predictions for W=2.5 GeV at five values of 
Q 2 = 1.6, 2.1, 2.6, 3.8 and 5 GeV 2 , top to bottom. 



section dependence on the angle $ between the electron 
and hadron scattering planes. The value of the <jlt is 
consistent with zero and the assumption that SCHC is 
valid for <f) production in this kinematic regime. How- 
ever, small deviations are still possible as shown by more 
accurate measurements of these parameters at HERA en- 
ergies [25( | . 

The ratio R = ul/o't has been determined from two 
projections of the angular decay distribution of the K + 
in the (/>-meson rest frame and under the hypothesis of 
SCHC. The measurement of r{$ gives R = 1.05±0.38 and 
the measurement of r\_ l gives a value of R = 0.72±0.30, 
the weighted average being R = 0.85 ± 0.24. This mea- 
surement can be compared to the value of R=1.25 pre- 
dicted by the JML 2-gluon exchange model. We note that 
these extractions, at least from Tqq, are relatively insen- 
sitive to the assumption of SCHC as shown in Ref. [25j. 

The measurements of R from this analysis and other 
world data are plotted as a function of Q 2 in Fig. [19] 3 . 
The data show that the ratio R is increasing as a func- 
tion of Q 2 , but <tl is still not dominant at these kine- 
matics. Using our measurement of R, we can compute 
the average longitudinal cross section for our data. The 
average cross section is given by u(Q 2 = 2.21GeV 2 ) — 
6.9 ±1.7 nb, which yields a longitudinal cross section 
a L = 4.5 ±1.1 nb. 



XI. SUMMARY 

</>-meson electroproduction was examined in the kine- 
matical regime 1.4 < Q 2 < 3.8 GeV 2 , 0.0 < t' < 3.6 



3 The W-dependence of R has been studied at HERA , which 
covers a very large range in W. 
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FIG. 19: (Color online) R = ol/ot vs. Q 2 for our data (solid 
circles), previous CLAS result (open circle), HERMES results 
(triangles) Cornell data (stars), ZEUS data (open diamonds) 
and HERA data (squares). The two determinations from the 
present analysis are separated for ease of viewing about the 
actual Q 2 value of 2.21 GeV 2 . 



gests that we have not yet reached the asymptotic regime 
where QCD factorization can be applied without sub- 
stantial corrections. 

The cross sections have a weak dependence on —t, 
which indicates that at this Q 2 , the photons couple to 
configurations of substantially smaller size than the tar- 
get. Our data provide a very precise measurement of the 
exponential slope at small cAr ~ 0.5 fm, which shows 
that we are probing very small distances, approaching 
about one third the size of the proton itself. A natu- 
ral explanation is that (f> production is dominated by the 
scattering of small size ss virtual pairs off the target. 
This conclusion is supported by the good agreement be- 
tween our data and the extension of the JML model from 
the real photon point (where it has been calibrated) to 
the virtual photon sector. It describes the interaction 
between this ss pair and the nucleon by the exchange of 
two dressed gluons. We conclude that these constituent 
degrees of freedom are appropriate for the description of 
</>-meson production at low W and Q 2 ~ 2-3 GeV 2 . 
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